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Abstract 

The use of multi-dimensional finite volume numerical 
techniques with finite thickness models for calculating 
aeroheating rates from measured global surface 
temperatures on hypersonic wind tunnel models was 
investigated. Both direct and inverse finite volume 
techniques were investigated and compared with the 
one-dimensional semi -infinite technique. Global 
transient surface temperatures were measured using an 
infrared thermographic technique on a 0.333-scale 
model of the Hyper-X forebody in the Langley 
Research Center 20 -Inch Mach 6 Air tunnel. In these 
tests the effectiveness of vortices generated via gas 
injection for initiating hypersonic transition on the 
Hyper-X forebody were investigated. An array of 
streamwise orientated heating striations were generated 
and visualized downstream of the gas injection sites. In 
regions without significant spatial temperature 
gradients, one-dimensional techniques provided 
accurate aeroheating rates. In regions with sharp 
temperature gradients due to the striation patterns two- 
dimensional heat transfer techniques were necessary to 
obtain accurate heating rates. The use of the one- 
dimensional technique resulted in differences of ±20% 
in the calculated heating rates because it did not 
account for lateral heat conduction in the model. 

Nomenclature 

c specific heat, J/kg/K 
Fo Fourier number 

h enthalpy-based heat transfer coefficient, kg/nT/s 
h re f reference heat transfer coefficient, kg/nr/s 


H enthalpy, J/kg 
k thermal conductivity, W/m/K 
L length, m 

1^ length of 1 st Macor insert 
Ly width of 1 st Macor insert 
1^ thickness of Macor inserts 
q" heat flux, W/m 1 2 
T temperature, K 
T 0 initial (pre-run) temperature, K 
t run time (measured from initiation of model 
injection), s 

ti time model is at test section centerline, s 
t C on correction time, s 
t e ff effective time, s 
S sum of square of errors 
x local coordinate along model length 
y local coordinate along model width 
z local coordinate along model thickness 

a thermal diffusivity, m 2 /s 
(3 thermal product, •Jpck , J/m 2 /K/s 0 ' 5 
p density, kg/m 2 
Subscripts 

aw adiabatic wall condition 
tl reservoir conditions 

t2 stagnation conditions behind normal shock 
w wall condition 
°° freestream static condition 

Acronyms 
SI semi -infinite 
FV finite volume 


Aerospace Technologist, Metals and Thermal Structures Branch, MS 396, Senior Member AIAA 
+ Aerospace Technologist, Aero thermodynamics Branch, MS 408A, 

: Aerospace Technologist, Aerothermodynamics Branch, MS 408A 
s Research Scientist, Aerothermodynamics Branch, MS 408A, Member AIAA 

This material is declared a work of the U.S. Government and is not subject to copyright protection in the United States. 


1 

American Institute of Aeronautics and Astronautics 



Introduction 

One of the standard techniques for calculating 
aeroheating rates from measured transient surface 
temperatures on hypersonic wind tunnel models has 
been the one-dimensional semi- infinite solid 
conduction method. This method was originally 
developed for use with discrete temperature sensors, 
such as thin film resistance gages. 1 The semi-infinite 
analysis assumes that heat doesn’t noticeably diffuse to 
the back surface of the model, so that the back surface 
temperature doesn’t vary. It is also assumed that there 
is no lateral heat conduction so that heat is conducted 
only in the direction of the model thickness from the 
top surface (surface exposed to aeroheating) to the back 
surface. This assumption is necessary when using 
discrete temperature sensors, since due to limitations on 
the number and spacing of the sensors that can be 
installed on the model multi-dimensional conduction 
cannot be considered. There are several limitations to 
the use of the 1-D semi-infinite method. The semi- 
infinite assumption is only valid for a specific time 
duration that is a function of model thickness and 
thermophysical properties. The solution assumes 
constant thermophysical properties. Furthermore, this 
technique is not suitable for regions with sharp 
temperature gradients across the surface, because lateral 
heat conduction may be significant compared to 
through -the-thickness heat conduction. 

To overcome limitations on the number of discrete 
sensors that can be installed on models, optical 
techniques such as infrared and phosphor thermal 
imaging have been utilized for providing global 
aeroheating data on hypersonic wind tunnel models. 

The phosphor thermography technique is routinely used 
at Langley Research Center for hypersonic wind tunnel 
aeroheating measurements. 2 The accuracy of infrared 
thermal imaging for transient surface temperature 
measurements on a flat plate model in subsonic flow 
has been evaluated by comparison with surface 
mounted thermocouples. 3 This technique has been used 
for aeroheating calculations on hypersonic wind tunnel 
models, 4 ' 5 and for in -flight surface temperature 
measurements on the Shuttle Orbiter. 6 ' 

The 1-D semi -infinite technique is the standard 
method used with thermal imaging techniques, even 
though the availability of whole -field surface 
temperature distributions lends itself to using multi- 
dimensional thermal analysis. The 1-D semi- infinite 
method yields a simple analytical solution, 8 while 
multi-dimensional solid conduction analysis requires 
numerical thermal analysis software and complete 
modeling of the geometry of the wind tunnel model. 
Nevertheless, the 1-D semi -infinite model has produced 
satisfactory results on blunt body configurations. 2 


Future vehicles, such as air-breathers, are slender with 
thin forebody/control surface leading edges. In 
addition, there is a critical need for more detailed 
studies of micro aerothermodynamic phenomena (such 
as vortices and protuberances). The 1-D semi -infinite 
technique is not adequate to accurately assess these 
issues with future vehicles. The present study is an 
attempt to investigate use of numerical multi- 
dimensional heat conduction models with global 
thermal data as the next level of data reduction 
methods. 

Numerical techniques have been previously used 
with discrete sensors to overcome some of the 1 -D 
semi -infinite limitations. The 1-D finite volume 
numerical technique has been used for calculating 
aeroheating rates from discrete sensors on hypersonic 
wind tunnel models. 9, 10 This technique numerically 
solves the governing diffusion equation using the 
measured temperatures on the top surface and assuming 
an adiabatic back face boundary condition. The main 
advantage of this technique is that it can use 
temperature -varying properties and it does not have the 
time limitation validity issues associated with the semi- 
infinite technique. An inverse 2-D finite element 
numerical technique has been used in conjunction with 
temperature measurements from a linear array of 
closely spaced thin film gages to obtain heating rate 
distributions in the region of shock impingement on a 
cylinder in Mach 6 How. 1 1 This technique resulted in 
higher peak heating rates at shock impingement 
locations, but the accuracy of the technique could not 
be determined due to limitations on the number of 
sensors used. 

In order to investigate multi-dimensional finite 
volume numerical techniques with global thermal 
imaging data, a series of runs were conducted as a 
subset of an overall experimental program for 
evaluating gas injection for tripping the boundary layer 
on a 0.333-scale Hyper- X forebody model. This 
tripping mechanism produces an array of streamwise 
orientated heating striations downstream of the gas 
injection sites. Various solid trip configurations had 
previously been investigated on this 0.333- scale 
Hyper-X forebody model. 12 The overall test program 
was a comprehensive study to investigate the 
effectiveness of gas injection for hypersonic boundary 
layer transition and included more than 80 wind tunnel 
runs. It used phosphor thermography for aeroheating 
measurements and included oil flow visualization runs. 
The runs reported here were a small subset of the 
overall test program and concentrated on evaluation of 
multi-dimensional solid conduction numerical methods 
in conjunction with infrared thermographic techniques 
to study the striation heating downstream of the gas 
injection sites. 
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The main objective of this investigation was to 
compare 1-D and multi-dimensional finite volu me 
methods for predicting heating rates, especially in 
regions with large spatial temperature gradients. Both 
direct and inverse finite volume techniques were 
investigated. 

Experimental Methodology 

Test Facility 

Aeroheating data were obtained on a 0.333-scale 
Hyper-X forebody model in the 20-Inch Mach 6 Air 
Tunnel of the LaRC Aerothermodynamics 
Laboratory. 13 This tunnel is a conventional blowdown 
tunnel with a 0.521 m x 0.508 m (20.5 x 20 in.) test 
section, and uses a two-dimensional contoured nozzle 
to provide a nominal freestream Mach number of 6. 

The tunnel uses dried, heated, filtered air as the test gas. 
A detailed description of the facility and its associated 
instrumentation is found in Ref. 13. Typical operating 
conditions for the tunnel are stagnation pressures of 
0.21-3.45 MPa, stagnation temperatures of 483-533 K, 
and free stream unit Reynolds numbers of 1.6 x 10 6 /m- 
25.6 x 10 6 /m. A bottom- mounted model injection 
system inserts models from a sheltered position to the 
tunnel centerline in approximately 1-2 s. Typical 
duration of aeroheating tests in this facility is 5- 10 s. 

Model 

A photograph of the 0.333-scale Hyper-X forebody 
model is shown in Fig. 1. This model is described in 
detail in Ref. 12, and was previously used for 
comparing the effectiveness of various hypersonic 
boundary -layer trip devices. 12 The windward forebody 
model consisted of three flat ramps that provided a 
series of discrete nonisentropic flow compression 
surfaces simulating the engine external inlet ramps. 

The model was tested at 2° angle of attack, and the first, 
second and third ramps provided additional 
compressions of 2.5°, 5.5°, and 3°, respectively. 

The forebody model, 71.12-cm long, was built 
from aluminum and had a detachable stainless-steel 
leading edge (nose radius of 0.254 mm) and 
interchangeable measurement surface inserts. The 
length of the leading edge was 12.7-cm, while the trip 
station was 18.84-cm from the model leading edge. 

The trip insert plates consisted of several stainless steel 
inserts with various orifice shapes and spacing for gas 
injection into the flow in order to trip the boundary 
layer. The specific trip insert plate used for this report 
consisted of a single row of 17 orifices, 0.25-mm in 
diameter with a 3. 18- mm spacing. The re maining flat 
ramp sections were designed to accommodate either a 


Macor or aluminum set of inserts. The Macor and 
aluminum inserts were used for heat transfer and flow 
visualization studies, respectively. The Macor inserts 
were 6.35-mm thick. The leading edge of the first 
insert plate, 10.87-cm long and 12.9-cm wide, was 
20.74-cm downstream of the model leading edge. 

Since phosphor thermography was the main 
measurement system used for obtaining heating rates 
for the overall experimental program, the Macor 
substrates were coated with a 0.025- mm thick mixture 
of phosphors suspended in a silica-based colloidal 
binder. A photograph of part of the model showing the 
leading edge, the trip insert plate with 17 orifices for 
gas injection, and the first Macor insert is shown in Fig. 
2 . 

Test Conditions 

The nominal reservoir stagnation pressure and 
temperature for the two runs reported in the present 
study are presented in Table 1. The reported flow 
conditions were based on the measured reservoir 
stagnation pressures and temperatures and recent 
unpublished tunnel calibrations. Both runs were 
conducted with the model at an angle of attack of 2°. 
Run 70 was a baseline test without gas injection into the 
boundary layer, while Run 71 consisted of gas injection 
into the boundary layer. The total pressure of the gas 
injected into the boundary layer was measured in a 
manifold directly below the orifices, and was 31.2 kPa. 
The total pressure was used to control the penetration 
depth of the jets compared to the boundary layer 
thickness. 

Test Technique 

Model surface temperatures were measured using 
an infrared imaging system with an un-cooled 
microbolometer-based focal plane array detector with 
320 x 240 detector elements. The imager has a field of 
view of 24° x 1 8°, and is sensitive to infrared radiation 
emitted in the 7-12 micrometer spectral range of the 
electromagnetic spectrum. An eight-inch diameter zinc 
selenide window with proper anti-reflection coating to 
provide transmittance of 0.98 in the 7-14 micrometer 
spectral range was used for optical access to the test 
section. The measured infrared radiation is converted 
to actual surface temperatures using the emittance of 
target. Ceramic models have high emittance values and 
do not require any special coating for infrared thermal 
measurements. But the Macor inserts used in this test 
had already been coated with phosphor coating to 
obtain global thermographic phosphor data in 
subsequent tunnel runs. The room-temperature, total 
(wavelength averaged), near normal emittance of the 


a machinable glass ceramic which is a registered 
trademark of Corning Incorporated 
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Macor inserts with the applied surface phosphor coating 
was determined to be 0.88 using a heated cavity 
reflectometer. 14 Joints between the Macor inserts 
corresponding to the location of ramp angle changes 
(Fig. 1) were used to geometrically locate the model in 
image mappings. The imager produces 30 frames of 
images per second, while the data acquisition hardware 
used was only capable of real-time digital storage of 
approximately 6 frames per second. Surface 
temperatures were measured on parts of the model 
leading edge, the entire first Macor insert, and part of 
the second Macor insert. But only the data from the 
first Macor insert was used in this study to compare 
various aeroheating calculation techniques. Model 
surface streamlines were obtained using the oil flow 
technique. Metal inserts were used on the model and 
were spray painted with black paint to enhance contrast 
with the white -pigmented oil. Surface streamlines were 
recorded using a conventional video camera. 12 


Solution Techniques 


1-D Semi -Infinite Method 

The standard 1-D semi-infinite solid conduction 
formulation with specified surface temperature 
boundary condition cannot be used with global optical 
techniques in hypersonic wind tunnels with model 
injection since detailed surface temperature data is not 
available during the time the model is being injected 
into the test section. Imagers are set such that they are 
focused on the model after the model has reached the 
centerline position in the test section. Furthermore, 
with optical techniques requiring ultraviolet excitation 
of the target, the model might not be completely and 
uniformly illuminated till it reaches the test section 
centerline. Therefore, the semi-infinite solid 
conduction formulation with convective boundary 
condition has been utilized for global thermographic 
measurement techniques. 2 The governing equation is 
the constant property 1-D heat conduction equation, 
with an enthalpy-based convective boundary condition 


-k—{0,t) = h[H aw -H w ] (1) 

oz 

H aw is the adiabatic wall enthalpy assumed to be equal 
to the total enthalpy (Hti in Table 1), and H w is the wall 
enthalpy evaluated at wall temperature [T w = T(0,t)]. 

The analytical solution is 


T(z, t) = T 0 + 
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where erfc is the complimentary error function. At the 
surface ( z=0 ), the solution reduces to 
T(0,t)-T o 
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The measured surface temperature is used in Eq.(4) to 
solve for A, which then yields the heat transfer 
coefficient, h, from Eq. (3.b). Since the formulation is 
for constant temperature thermophysical properties, 
calculations are made using temperature -averaged 
properties. 2 The above calculations are performed for 
every point (pixel) in the image. All that is needed to 
compute heat transfer coefficient distributions over the 
model is the initial surface temperature distribution, one 
test image, and the corresponding time. This technique 
assumes that the heat transfer coefficient distribution is 
time invariant. This assumption is not valid when the 
model is traversing through the test section’s boundary 
layer. In order to overcome this difficulty, the injection 
process is modeled as a delayed step heating, with the 
delay in the step corresponding to the time when the 
model has traversed halfway through the tunnel 
boundary layer.” The time it takes from the initiation of 
the injection process until the model traverses halfway 
through the tunnel boundary layer edge is referred to as 
the “correction time”, t cor r, and is subtracted from the 
run time to obtain an effective time to be used in the 
above equations 2 

t e ff~t~tcorr (^) 


where t is the actual run time which is triggered when 
the injection system first starts to move. The 
calculation of t corr is described in detail in Ref. 2. This 
technique has been validated by aeroheating 
calculations on a 50.8 -mm diameter hemisphere in 
Mach 10 flow, 2 and has been successfully applied to 
typical blunt body configurations such as X-34 2 and X- 
33 15 with favorable comparisons with computational 
fluid dynamic (CFD) calculations. 


Finite Volume Numerical Technique 

As mentioned previously only the data from the 
first Macor insert are discussed in this study in order to 
simplify geometrical considerations and concentrate on 
the basic aeroheating calculation issues. Referring to 
Fig. 3 for local coordinate system and geometrical 
definitions, the governing 3-D diffusion equation and 
initial and boundary conditions used for the Macor 
insert were 
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These equations are valid for t > ti, with ti 
corresponding to run time when the model has reached 
the test section centerline, since detailed surface 
temperature data are only available after this time. For 
simplicity it was assumed that adiabatic boundary 
conditions existed at the edges of the Macor insert (y = 

0, y = Ly, x = 0, and x = L*), Eqs. (6.e)-(6.h). Actually, 
heat is transferred from the edges of the first Macor 
insert to the adjoining surfaces on the model (metallic 
surfaces or second Macor insert), but the exact 
modeling of these boundary conditions was not 
necessary, since the imposed boundary conditions do 
not significantly affect heat transfer in the region of 
interest away from the boundaries. It was assumed that 
the model back face (z = LJ is adiabatic [Eq. (6.d)], 
while the top surface has a specified transient 
temperature distribution, f(x,y,t), measured using the 
optical thermographic technique [Eq. (6.c)]. The Macor 
inserts were placed on the model using narrow spacers, 
resulting in air gaps between the back face of the inserts 
and the top of the metallic model. Therefore, the 
adiabatic back face boundary condition was utilized. 

This assumption is certainly valid during the time the 
semi -infinite assumption is valid. Furthermore, after 
the back face temperature begins to rise, the heat losses 
from the back face are typically low enough that the 
assumption of an adiabatic condition may not introduce 
large errors The initial temperature distribution through 
the thickness of the model at time t = t ; is represented 
as g(x,y,z). This initial temperature distribution was 
not known for z > 0, and the z variation was 
approximated at each pixel location by the analytical 
solution of 1-D semi-infinite conduction model, Eq. (2), 
at time ti. Any errors in this initial temperature 
distribution will be corrected after a few time steps of 
the numerical solution due to the diffusive nature of 
heat conduction. 


The governing equations for the 2-D and 1-D finite 
volume techniques are subsets of the set of equations 
provided in Eqs. (6.a>(6.h). For the 1 -D case the 
dependence on x and y coordinates and the 
corresponding boundary conditions in these directions 
are eliminated. For the 2 -D case the dependence on x 
coordinate and the corresponding boundary conditions 
in this direction are eliminated. The solution to this set 
of equations constitutes the direct solution method. The 
equations were solved using a finite volume numerical 
method with variable thermophysical properties. The 
transient problem was solved using unconditionally 
stable implicit time marching techniques: 16 Crank 
Nicholson method for the 1-D case and Alternating- 
Direction Implicit (ADI) method for the 2-D and 3-D 
cases. The governing equations and the derivative 
boundary conditions were formulated such that second- 
order accuracy could be achieved with truncation errors 
on the order of At 2 and Ax 2 , where At and Ax are the 
time step and spatial nodal spacing, respectively. For 
the results reported here 21 nodes were used in the z- 
direction, and 128 by 159 nodes in the y and x 
directions, respectively. The latter was dictated by the 
number of infrared image pixels across the top surface 
of the first Macor insert. Solutions obtained using 5 1 
and 21 nodes in the z direction were similar, thus 
verifying that the nodal spacing used was sufficient. 

The time step used was dictated by the digital storage 
interval of infrared images, which was approximately 
0.17 s. Upon numerical solution of the governing 
equations, the convective heat transfer coefficient was 
calculated at each instant of time from 
3 T 

-k — ( x, y,0, t) = h (x ,y,t ) 

& (6.i) 

[H aw -H w (x,y,t)] 

For the inverse method, the conservation of energy 
equation, Eq. (6. a) was solved subject to all the initial 
and boundary conditions stated above, except that the 
specified temperature boundary condition on top 
surface, Eq. (6.c), was replaced by the convective 
boundary condition in Eq. (6.i), and it was assumed that 
the convective heat transfer coefficient was independent 
of time. Since the convective heat transfer coefficients 
were unknown, the Levenberg-Marquardt method for 
nonlinear parameter estimation 1 7 was utilized to 
estimate their values. The inverse technique is an 
iterative numerical scheme where a matrix h of 
estimates representing h(x,y) is assumed; the numerical 
solution to the equation using this boundary condition 
and previously described boundary and initial 
conditions are marched in time, and surface temperature 
distributions are calculated. The comparison of 
measured and calculated transient surface temperatures 
and the Jacobian matrix representing the derivatives of 
surface temperatures with respect to the heat transfer 
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coefficients are used to develop better estimates of h. 
This iterative procedure is based on the minimization of 
the ordinary least squares norm given by 

S = 'L[T m (x , y,0,t) -T p (x, y, 0, t) ] 2 (7) 

where T m and T p are measured and calculated surface 
temperatures, respectively. The 1-D inverse method 
required estimating one heat transfer coefficient at each 
measurement point on the image. The 2-D and 3-D 
inverse methods required simultaneous estimation of 
either an array of heat transfer coefficient along a 
specific axis, or a matrix of heat transfer coefficients 
over the entire image, respectively. The computational 
effort increases tremendously as analysis is extended 
from 2-D to 3-D. 

For the inverse method, transient solution could be 
obtained at time steps smalle r than the time increment 
at which surface temperature distributions were 
available, but the comparison of predicted and 
measured temperatures, required in Eq. (7), was 
performed at times when surface temperature data were 
available. 

Results and Discussion 

Results for the baseline run without gas injection 
(Run 70) are presented to compare various techniques 
for calculating aeroheating rates. Then the results for 
the run with gas injection (Run 71) are presented. 

Macor thermophysical properties were obtained from 
published results in Ref. 9. The convective heat 
transferrates were non-dimensionalized with respect to 
a reference heat transfer coefficient, h re f, corresponding 
to the Fay and Riddell 1 8 stagnation-point heating on a 
101. 6- mm radius sphere. The standard data reduction 
procedure for conducting aeroheating tests using global 
thermographic techniques was used. The model was 
injected into the test section just prior to establishing 
tunnel flow and a “pre -run” image of initial surface 
temperature distribution across the model was obtained. 
Then the model was retracted. After the flow was 
established, the model was injected into the test section, 
and the start of the injection process triggered the 
infrared imaging digital recording process at 0.17 s 
intervals. Only the surface temperature and aeroheating 
results on the first Macor insert are presented. 

Run without Gas Injection 

The run without gas injection provided uniform 
surface temperature and aeroheating distributions 
across most of the surface area of the 1 st Macor insert. 
Therefore, this run was utilized to assess the accuracy 
of 1 -D and multi-dimensional direct and inverse finite 
volume methods. The initial (pre -run) average surface 
temperature on the first Macor insert was 303.7 K with 
a standard deviation of +0.28 K. Run time of 1.37 s 


corresponded to the time when the first image was 
digitally recorded after the model had reached the test 
section centerline. The spatial distribution of surface 
temperature rise compared to initial temperature along 
the centerline of the forebody model (y = Ly/2) at 
various times throughout the run is shown in Fig. 4. 
Data are plotted versus the non-dimensional length, 
x/Lx. The surface temperatures rose only 6.1 K after 8.5 
s for this low aeroheating condition; i.e., a slender body 
at a low angle of attack (2°) and at unit Reynolds 
number of 7.3 x 10 6 /m. The adverse pressure gradient 
provided by the second ramp induced a local flow 
separation near the end of the first Macor insert, as is 
evident from the sharp decrease of surface 
temperatures, especially at later times. 

The spatial variation of non-dimensional 
convective heat transfer rates along the model 
centerline (y = Ly/2) calculated using various 1-D 
techniques is shown in Fig. 5: semi-infinite, direct and 
inverse finite volume results. The semi-infinite and 
direct finite volume results shown are at t = 2.58 s. 

Semi -infinite results are presented both with and 
without time corrections, Eq. (5). Originally, the semi - 
infinite results were calculated using the standard 
technique 2 with the appropriate time correction, which 
for the runs reported here was 0.89 s. But comparison 
of the results with other experimental results showed 
that these results were over-predicting the other results 
significantly. Therefore, semi-infinite results were also 
calculated with t corr =0 s which better agreed with other 
results. It appears that the use of time correction with 
semi -infinite techniques is not applicable to wind tunnel 
tests on slender bodies at low heating rates. Therefore, 
for the direct and inverse finite volume numerical 
methods, the initial temperature distribution at 
t=t i=1.37 s, was calculated from the semi-infinite 
solution, Eq. (2), with t corr =0. The inverse results 
shown in Fig. 5 were calculated using temperature data 
from 6 thermal images with 0.17 s time steps with the 
first image at 1.37 s. Using 6 thermal images with 0.34 
s time steps and 12 thermal images at 0.17 s time steps 
yielded similar results. Overall, there was good 
agreement between direct and inverse finite volume 
results. The semi -infinite results with no time 
correction slightly over-predicted the finite volume 
results, while the semi -infinite results with time 
correction were significantly higher than all the other 
results. 

The temporal variation of the heat transfer 
coefficients spatially averaged over the interval 0 <x< 

1^ along the model centerline (y=Ly/2) using semi- 
infinite technique with and without time correction and 
direct and inverse finite thickness numerical results are 
shown in Fig. 6. The semi -infinite technique with time 
correction over-predicted the heating rates significantly 
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at early times, then, asymptotically approached a 
constant value. The semi-infinite technique with no 
time correction slightly over-predicted the heating rates 
at early times, and quickly leveled off. The semi- 
infinite results are very sensitive to correction time, 
especially at early times. Even though the semi -infinite 
results with time correction provided correct results 
after t = 6s, the validity of the semi -infinite assumption 
at this time is questionable. The semi-infinite 
assumption for a convective boundary condition is valid 
for 19 

at 

Fo = — <0.097 (8) 

L \ 

With Macor thermophysical properties at room 
temperature and a model thickness of 6.35 mm, the 
semi -infinite assumption is only valid for t < 4.9 s. So, 
even though reducing data at later times (t > 6 s) might 
produce accurate results, but the basic assumption used 
in the analytical solution is violated. The direct finite 
volume results also slightly over-predicted the heating 
rates at early times, and then slightly oscillated around a 
constant value at later times. The inverse finite volume 
technique generated time independent results which all 
the other techniques’ results eventually approached. 

The direct and inverse 1-D finite volume results 
described previously were obtained by using the semi - 
infinite solution without time correction as the initial 
temperature distribution (t = 1 1 ) through the thickness of 
the model. The sensitivity of the inverse 1-D finite 
volume result to this initial temperature distribution was 
investigated by comparing this result to one when time 
correction was used for the initial temperature 
distribution at a point at the center of the 1 st Macor 
insert, x = L^/2, y = Ly/2. The resulting convective heat 
transfer coefficients with and without time corrections 
for the initial temperature distribution were 0.12 and 
0.1 1, respectively. Therefore, the inverse finite volume 
technique appears to be fairly insensitive to the initial 
temperature distribution used. 

The variation of convective heat transfer 
coefficients with y/Ly at x/L* = 0.75 calculated using 
one, two, and three dimensional direct finite volume 
formulations are presented in Fig. 7. The 2-D results 
were calculated at each longitudinal position across the 
width of the model from y/Ly = 0 till y/Ly =1. For the 
two-dimensional case the dependence on jt coordinate 
and the corresponding boundary conditions in this 
direction were eliminated. The three sets of results 
coincided with each other except at the edges (y/Ly = 0 
and 1), where the 2-D and 3-D results slightly varied 
from 1-D results due to the imposition of the adiabatic 
boundary conditions at the edges for the multi- 
dimensional formulations. This showed that the 
baseline case had no in -plane gradients of a magnitude 


affecting the calculations, and that for this case 1-D 
analysis was sufficient. 

Run with Gas Injection 

A run was conducted with gas injection into the 
boundary layer for tripping the boundary layer flow. A 
photograph of the oil flow visualization for a run with 
similar flow conditions is shown in Fig. 8. The orifices 
for gas injection, and the streamlines corresponding to 
18 vortices downstream of the gas injection sites are 
visible in the photograph. Each gas injection orifice 
produced two vortices, with the two vortices between 
adjacent orifices merging together, thus, producing 18 
vortices for this 17 -orifice configuration. A contour 
plot of the corresponding non-dimensionalized 
convective heat transfer coefficients calculated using 
the inverse 1-D finite volume technique is shown in 
Fig. 9. The traces of 18 streamwise oriented heating 
striations downstream of the gas injection sites are 
visible in the contour plot. The heating results were 
consistent with oil flow results in a qualitative sense. 

It was believed that in the vicinity of the vortices 
with high spatial temperature gradients 1-D heating 
results might not be accurate, therefore, multi- 
dimensional finite volume techniques were utilized. 
First, direct 2-D and 3-D direct finite volume 
techniques were compared to determine which 
formulation would be more appropriate for the given 
run conditions. The 2-D results were calculated at each 
longitudinal position across the width of the model 
from y/Ly = 0 till y/L, =1 (the dependence onr 
coordinate and the corresponding boundary conditions 
in this direction were eliminated). The comparison of 
the results at the longitudinal position of x/L x = 0.25 is 
shown in Fig. 10. The heating peaks related to the 18 
vortices are visible between y/Ly of 0.2 and 0.8. The 
heating rates at some vortex locations are as high as 3 
times the heating rates in the undisturbed regions. 
Furthermore, it is seen that there is no difference 
between the heating results predicted using either the 2- 
D or 3-D analysis, with the root mean square difference 
between the two sets of results being 1%. Therefore, it 
was determined that a 2-D analysis is sufficient for 
obtaining heating rates for this run. It was also 
desirable to use the inverse method since it yields time- 
invariant heat transfer coefficients. Therefore, 2-D 
inverse finite volume analysis was used for the rest of 
the study. The comparison of 1 -D and 2-D inverse 
finite volume results for the longitudinal positions of 
x/LxOf 0.25 and 0.5 are shown in Figs, ll.a and ll.b, 
respectively. The 1 -D results under-predict the actual 
heating rates at the peak heating locations, and over- 
predict the heating rates in the valleys between two 
adjacent peaks. Lateral heat conduction smears the 
surface temperature distribution in regions with sharp 
heating gradients; heat gets transferred from the peaks 
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to the valleys, causing reduction of temperature at the 
peak and increase of temperature at the valleys. When 
1-D heat transfer is utilized with the measured 
temperature distributions, the predicted heating rates 
are lower at the peaks and higher at the valleys. Only 
multidimensional, in this case 2-D, heat transfer 
analysis can predict the correct heating magnitudes, 
since it accounts for lateral heat conduction. Any noise 
in the measured surface temperatures gets amplified in 
the estimation of the heating rate, and regularization 
techniques have previously 11 ' 20 been used utilized with 
inverse techniques to overcome this problem. In the 
present study the noise in the measured surface 
temperatures did not pose a major problem for heating 
rate calculations, and no regularization technique was 
used. The difference between heating rates calculated 
using 1-D and 2-D methods along the y axis at 
longitudinal positions of x/Lx of 0.25 and 0.5 are shown 
in Figs. 12. a and 12.b, respectively. The 1-D method 
under-predicts and over-predicts the heating rates at the 
peaks and valley by almost 20%. 


Concluding Remarks 

Global transient surface temperatures were 
measured using an infrared thermographic technique on 
a 0.333 -scale model of the Hyper-X forebody in the 
Langley Research Center 20 -Inch Mach 6 Air tunnel. 

In these tests the effectiveness of vortices generated via 
gas injection for initiating hypersonic transition on the 
Hyper-X forebody were investigated. An array of 
streamwise orientated heating striations were generated 
and visualized downstream of the gas injection sites. 

The striation patterns were attributed to the 
development of embedded boundary layer vortices. 

The use of multi-dimensional numerical finite volume 
numerical techniques for calculating aeroheating rates 
from measured global surface temperatures was 
investigated. Both direct and inverse finite volume 
techniques were investigated and compared with the 1 - 
D semi -infinite method. In regions without significant 
spatial temperature gradients, 1-D techniques provided 
accurate aeroheating rates. The use of time correction 
with 1-D semi -infinite results did not produce 
satisfactory results for this low heating test on a slender 
body. The 1-D inverse finite volume technique resulted 
in time -invariant heating rates which all the other 
techniques’ results eventually approached. In regions 
with sharp temperature gradients due to the striation 
patterns 2-D finite volume techniques were necessary to 
obtain accurate heating rates. The use of the one- 
dimensional technique resulted in differences of ±20% 
in the calculated heating rates in the vicinity of vortices 
because it did not account for lateral heat conduction in 
the model. 
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Figure 2. Photograph of trip insert plate with 17 orifices 
and first Macor insert on the Hyper-X forebody model. 
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Figure 3. Coordinate system and geometrical 
definitions for the first Macor insert. 



Figure 4. Spatial variation of temperature along first 
Macor insert centerline (y= Ly/2) at various times. 



Figure 5. Spatial variation of aeroheating rate along 
first Macor insert centerline using various techniques. 



t, s 


Figure 6. Temporal variation of spatially averaged 
aeroheating rate (0<x/L,< 1 -0) along 1 st Macor insert 
centerline using various techniques. 



Figure 7. Comparison of 1 -D, 2-D and 3-D direct finite 
volume techniques at x/L x =0.75 at t= 2.58 s for baseline 
test without gas injection. 



Figure 8. Oil flow photograph 
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h/h ref 

0.4 

0.385 
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0.355 

0.34 
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0.31 

0.295 

0.28 

0.265 

0.25 

0.235 

0.22 

0.205 

0.19 

0.175 

0.16 

0.145 

0.13 

0.115 

0.1 



b) 


Fig 9. Spatial variation of aeroheating rates for test with 
gas injection using 1-D inverse finite volume technique. 


Figure 11. Comparison of 1-D and 2-D inverse finite 
volume aeroheating results for test with gas injection at 
x/L x = a) 0.25, b) 0.5 



y/L y 


Figure 10. Comparison of 2-D and 3-D direct finite 
volume results at x/L*=0.25 for test with gas injection 



y/L y 

a) 



a) 



y/L y 

b) 


Figure 12. Relative difference between 1-D and 2-D 
inverse finite volume aeroheating results for test with 
gas injection at x/I^= a) 0.25, b) 0.5 
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